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Abstract 

We develop a model for spiral galaxies based on a nonlinear realiza- 
tion of the Newtonian dynamics starting from the momentum and mass 
conservations in the phase space. The radial solution exhibits a rotation 
curve in qualitative accordance with the observational data. 
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1 Introduction 

Galactic rotation curves are one of the most important evidences favoring the 
dark matter scenario in astrophysics. Actually dark matter was proposed by 
Zwicky 1 in order to accommodate observed velocities in galaxies and galaxy 
clusters. Since then others explanations was proposed as modifications of New- 
tonian gravity [2] or modifications in General Relativity |3J . All these approaches 
have in common, from the modeling point of view, just one element: non- 
linearity. So, is natural asks the question if non-linear effects of some sort are 
responsible for the discrepancy observed in galaxies and galaxy cluster dynamics. 
Recently, in the context of General Relativity, some authors have investigated 
this point [4]. 

Here, we will propose an alternative paradigm where Newtonian gravity is 
maintained on its solid basements but a non-linear model of a galaxy is build. 
We solve the non-linear equations in two simplified cases and calculate the 
resulting galaxy rotation curves, showing qualitatively that flat curves can be 
obtained in a given region of the parameter space. 
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2 The General Model 

We focus on spiral galaxies described as a fluid with mass distribution in the 
phase space, (x, v, t) which is related to the ordinary mass distribution by 



p (x, t) — J dv * (x, v, t) 



where m* the mass for a typical star, the dynamics is governed by the Poisson 
equation coming from the second Newton law, 

V 2 V(x,t) =47rGp(x,i) 

but restricted to the mass conservation, 

_ () ^ a* a* , 9* _ 

eft <9x 9v 

Using a Hamiltonian description, x = v, v = — VV, we obtain, 

— + v-v*-vy- — = o 

at av 
V 2 V (x,t) = 4nG y dv*(x,v,t) 

In general, this is a set of intcgro-diffcrential non-linear coupled equations. 
2.1 Axial Symmetry 

Let us restrict the model only to the disc of the galaxy in the static regime. 
So, using the axial symmetry of this system, the mass conservation equation 
become: 

^dH> + ^ + .d^ fdVd^ + dVd^ + dV dV\ _ Q 

dr d(j) dz \ dr dr d<j> d<j> dz dz J 

This equation can be solved using the Method of Characteristics. Assuming 
a constant angular velocity, this is equivalent to a set of ordinary differential 
equations, 

dr dr dz dz d<p dty 

~T = dV = ~T = dv ' ~T = ' ~aT = 
r — gjr z <f> dt 

whose uniparametric family of solutions is 

* = , E = l - (r 2 + r 2 ^ 2 + z 2 ) + V (r, z) 

Experimental data shows that the mass distribution is Gaussian in the ob- 
served velocities, therefore is natural to choose the mass distribution W to be a 
Boltzmann distribution in the total energy: 



# (E) = ^ o e~ 0E = *o exp ( -/3 



l - (r 2 + r 2 2 + i 2 ) + V (r, z) 
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The integration over the velocities can be done, in order to obtain the func- 
tional dependence of the mass distribution over the gravitational potential 

// n \ 3/2 
drd<pdzr^ (E) = * ( ~j J exp (-0V (r, z)) = p e~ mr ' z) 

It means that the Poisson equation now is a non-linear partial differential 
equation, given by 

or z r or oz z 

3 Variation in the Height 

Let us to take an over simplified case, assuming that the field vary only with 
the height to the plane of the disc, 

^ = AkG P {z) = ATtG Po e- fiv{ ^ 

Using an integration factor ^ and the boundary conditions V (0) = dV J^ = 
this equation can be directly integrated, 

V(z) = | In cosh (—') , zq = (27rGp /?r 1/2 



Po 



13 \z 
The density profile in this case is 



cosh 2 
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So far, it is a reasonable model, since the disc predicted is thin. 



4 Radial Variation 

Let us to take another over simplified model assuming variation only in the 
radial direction, 

dr 1 r dr 

Multiplying both sides by r and integrating, we find 

r^- — r a c = iirGpo [ fe~^ v ^dr 
dr J ro 

where r is the radius of the galaxy core, and c = dV ^°^ ■ Performing a second 
integration, we arrive in a Volterra second order integral equation: 

V (r) = V + r cln | — ) + 4ttG Po f dr f ^e'^^df 

VoJ Jro Jr r 
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To solve it, we apply the Piccard's method of successive approximations: 
( r ) = V<® (r) + AnGpo f df f V^^df 



V {0) (r) = V + r cln f 



Therefore, the first order iterative solution is 

2 



AirGpoe 



2 - (3r c ) \r ) 



5 Galaxy Rotation Curve 



Assuming a virial balance, ^- = ^ — > v = y^^r, the galaxy rotation curve in 
this model is U3 = (3r c , (r G )~ 2 = AnGp , A = ^ , a = £M 



In our illustrative example, rt = /3 = 1, we have the behavior illustrated in 
Fig. [I] 

So, if /3 > 2, 



2-/3 



which is constant. 



6 Final Remarks 

A non-linear Newtonian model for the galaxy disk was constructed by introduc- 
ing a solution of the mass conservation on the phase space. 

We illustrate by a direct example that the non-linear character of the gravi- 
tational field is a key feature to understand galaxy rotation curves, even in the 
Newtonian case. 
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Figure 1: Parametric rotation curve. Velocity is given in units of 1//3 1 / 2 and A 
is the distance to center in adimensional units. (3 is the parameter dictated by 
the boundary condition at the core. 



